Giant Nernst effect and bipolarity in the quasi-one-dimensional metal, Li .9Mo 6 Oi7 
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The Nernst coefficient for the quasi-one-dimensional metal, Lio.gMoeOi7, is found to be among 
the largest known for metals (y ~ 500 pV/KT at T ~ 20 K), and is enhanced in a broad range 
of temperature by orders of magnitude over the value expected from Boltzmann theory for carrier 
diffusion. A comparatively small Seebeck coefficient implies that Lio.gMosOi? is bipolar with large, 
partial Seebeck coefficients of opposite sign. A very large thermomagnetic figure of merit, ZT ~ 0.5, 
is found at high field in the range T ~ 35 — 50 K. 



The Nernst effect, a transverse electric field generated 
by a longitudinal temperature gradient in the presence 
of a perpendicular magnetic field, has been known since 
its discovery in the 19th century in the semimetal bis- 
muth [1]. Conductors with large Nernst coefficients are 
potential candidates for use in cryogenic refrigeration [2] , 
but in spite of several decades' research and a demand 
for such devices at temperatures below 100 K (e.g., for 
infrared detectors), none have yet proved practical. In 
recent years, Nernst-effect measurements have received 
renewed attention as a means for probing the pseudogap 
phase of cuprate superconductors [3, 4] and novel ground 
states of other correlated-electron systems [5] , prompting 
a reassessment of the physics leading to large Nernst co- 
efficients in metals. Emphasizing the contribution from 
diffusing charge carriers, experiment and theory agree [5] 
that the Nernst-coefHcient magnitude is set by the ratio 
of charge carrier mobility to Fermi energy. Here we re- 
port observations of very large Nernst coefficients in the 
quasi-one-dimensional (Q1D) metal, L10.9M06O17, rival- 
ing the largest known for bismuth and graphite. This 
material possesses neither a very large mobility nor a 
small Fermi energy. We establish that conduction along 
the Q1D chains of this material is bipolar, with very large 
partial thcrmopowers of opposite sign that predominate 
over conventional carrier diffusion in a broad tempera- 
ture range. The observations offer new insight into the 
origin of large thermoelectric coefficients in bulk conduc- 
tors and expand the criteria for viable materials in ther- 
momagnetic cooling applications. 

Lio.gMo60i7, a low-temperature superconductor (T c w 
2 K) first synthesized and studied in the 1980s [6, 7], 
has attracted interest more recently for its quasi-one di- 
mensionality and Luttingcr- liquid candidacy [8-11]. This 
compound is unusual among inorganic and organic Q1D 
compounds in its absence of a conventional density-wave 
transition [12, 13] (cither charge or spin) throughout a 
broad temperature range, T > T c . An upturn in its re- 
sistivity below Tm ~ 30 K may be associated with local- 
ization, dimensional crossover or the development of un- 
conventional (e.g., electronically-driven) charge density- 
wave order [12-14]. Though the T — > ground state (su- 



perconducting, metallic or insulating) is sensitive to heat 
treatment and can vary among as-prepared crystals, [15] 
the anisotropy and temperature variation of electrical 
properties for T > 5 K, the regime of study here, are 
robust and reproducible [12-17]. 

Single-crystal growth of Lio.9MoeOi7 is described in 
detail elsewhere [7, 13]. Crystals had typical dimensions 
0.05 — 0.15 mm along a and 0.5 — 2 mm in the b-c plane. 
Longitudinal and transverse voltages were measured with 
Au leads attached with silver epoxy. A 25 /xm-diameter 
chromel-constantan differential thermocouple monitored 
the temperature difference. The Nernst signal at each 
field was determined from the slope of linear Vh-AT x 
curves (with Vh the field-odd transverse voltage). Lat- 
tice constants and oxygen content for six crystals studied 
and additional details on the transport measurements are 
described in the Supplemental Material [18]. 

Complete sets of transport measurements (resistiv- 
ity, thermopowcr, thermal conductivity, Hall and Nernst 
coefficients) were measured on two single crystals of 
Lio.gMo60i7 (LMO), designated A2 and C, with electric 
or heat currents along the Q1D Mo-0 chains (crystallo- 
graphic b axis; inset Fig. 1). The Nernst coefficient of 
specimen C was also measured with heat current along 
the c axis. Thermal and thermoelectric measurements 
(heat flow along b) were also performed on a third crys- 
tal, designated D. For heat flow along the x direction 
(temperature gradient, V X T) and magnetic field along z, 
the Nernst coefficient is defined -E y /(V X TB Z ). 

With magnetic field along the least conducting a axis, 
Vcb and Vbc were determined for specimens A2 and C 
and v c b for specimen D as shown in Fig. 1(a) [18]. The 
observed ratio, iy c b/\^bc\ ^3 — 6, reflects the resistiv- 
ity anisotropy, p c / Pb, and agrees well with measurements 
on similar crystals [16, 17]. v(T) in nonmagnetic metals 
may have contributions from charge-carrier diffusion and 
from phonon drag, the latter arising from the interaction 
of the carriers with a non-equilibrium population of heat- 
carrying phonons. The diffusion term can be expressed as 
[5, 19], v/T ~ (ir 2 /3)(k B /e)(p/T F ) = 283(p/T F ), where 
p, is the carrier mobility and Tf the Fermi temperature. 
This expression approximates v/T quite well at low T 
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FIG. 1: (color online) (a) Nernst coefficient, v{T), for crystals A2 (dashed lines), C (circles and squares), and D (triangles). 
For crystal C: v c t and Vbc determined with primary heat current along b (solid circles) and c (solid squares), respectively. The 
specimen, originally with thickness 110 /im along a, was mechanically thinned to 55 /jm and remeasured with primary heat 
current along b (open circles) [18]. (b) u/T vs. T for crystal C, compared to bismuth [Ref.'s 20 (solid curve) and 21 (dashed 
curve)] and graphite (Ref. 22). The dash-dot-dot curve labeled "283/i/TF" represents the carrier-diffusion contribution from 
Boltzmann theory, determined from the carrier mobilities and estimated Fermi temperature [18]. (c) Figure of merit as a 
function of magnetic field for crystal C at several temperatures, (d) zero-field thermal conductivities along the b and c axes, 
inset: same data plotted as kT vs. T to emphasize the large constant term in Kb(T). Main Inset: Perspective view of the LMO 
monoclinic unit cell viewed along the b axis; Q1D Mol-0-Mo4 zigzag double chains composed of MoC>6 octahedra (highlighted 
bonds) are part of planar Mo-O networks whose two dimensionality is broken by the Li ions. 



(within a factor of 3 — 5) for Bi, graphite and a variety 
of other multi-band, correlated, and/or low-dimensional 
metals [5], over six orders of magnitude in fJ,/Tp [18]. 
The v/T data for LMO [Fig. 1(b)] exceed an estimate of 
this diffusion contribution (dash-dot-dot curve) by four 
orders of magnitude at low T, and remain significantly 
larger to room temperature. The expression is strictly 
valid for a linear energy dependence of the Hall angle at 
the Fermi level. Though strong energy-dependent scat- 
tering, as may arise from electron correlations or incipient 
density- wave order, or proximity to a quantum critical 
point [5], can enhance this energy dependence and hence 
v/T, it is very unlikely such effects can account for this 
very large discrepancy [19]. Wc conclude that a mech- 
anism other than carrier diffusion predominates in v{T) 
over the entire temperature range. 

Consider the potential of LMO for application. For 
operation as an Ettingshauscn cooler, with electric cur- 



rent applied along the b axis generating heat flow along 
c, the relevant adiabatic thcrmomagnctic figure of merit 
is [2], Z bc T = N bc N cb T / (p b n c ) where N = vB and k c 
is the c-axis thermal conductivity. Fig. 1(c) shows the 
field dependence of Z bc T along with n c (T) [Fig. 1(d)]; 
the latter is independent of field. We see that ZT at 
35 K reaches 0.5 for B=9T with no sign of saturation, 
among the highest known values for any material [23, 24]. 
Viable thermomagnetic materials for applications should 
have ZT approaching 1 at fields achievable using small 
permanent magnets (B < IT). It is conceivable that the 
inter-chain k c can be reduced (e.g. through the intro- 
duction of mass disorder via chemical substitution) to 
enhance ZT at lower field. 

The Seebeck coefficient or thermopower (TEP), S = 
—E x /\/ x T, reveals an important aspect of conduction 
in LMO. Crystals C and D with the largest v(T) have 
TEP values (Fig.'s 2 and 3) that are modest in magni- 
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FIG. 2: (color online) Seebeck coefficient (S), resistivity, 
and Hall coefficient for specimen C. Inset: Fermi surface for 
Lio.gMosOiT (adapted from Ref. 9). 
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FIG. 3: (color online) Zero-field &-axis Seebeck coefficient 
for four Lio.gMo60i7 specimens. Combined results of x-ray, 
compositional, and Hall measurements indicate that speci- 
mens with large peaks (Al, A2) are more oxygenated and 
particle-hole asymmetric (R h ^> R e and a h < a e ) compared 
to those with linear-T thermopowers (C,D,E) [18]. Inset: 
Hole and electron partial thermopowers computed using the 
equation [25], the measured Nernst coefficient [v c b), and two- 
band model parameters; error bars (25%) are propagated from 
those of the Hall data [18] . 



tude. Consequently, the ratio of transverse to longitu- 
dinal thermoelectric fields, E y /E x = vB/S, exceeds 10 2 
at 40 K; the electric field is essentially perpendicular to 
the temperature gradient. The implication is that LMO 
is compensated, with a cancellation of very large hole 
and electron Seebeck coefficients (opposite in sign). This 
can be seen more explicitly from the expression for v 
in a bipolar, anisotropic conductor [25] which includes a 
bipolar term involving the difference, S h — S e , of the hole 
and electron partial Seebeck coefficients along the heat 
flow. For perfect compensation, v c \, ~ (l/2)/ifc (S£ — S'f), 
while Sb = (l/2)(S b h + S' b : ). Bipolarity is plausible given 
the band structure [8, 9, 26] composed of two nearly- 
degenerate, approximately 1/2-filled, Q1D bands of Mo 
d xy parentage crossing the Fermi energy; there are two 
pairs of slightly warped Fermi surface (FS) sheets in the 
b*-c* planes (inset, Fig. 2). Each sheet, with positive 
and negative curvature, can contribute electron and hole 
character if the mean free path is anisotropic [27]. Bipo- 
larity is well known to enhance the Nernst effect and is 
a prerequisite for useful thermomagnetic materials. 

Supporting this bipolar picture, the TEP of other crys- 
tals with higher oxygen content, like Al and A2 (Fig. 3), 



exhibit very large, positive maxima at T = 50 — 70 K 
that augment the linear-T form (specimens C, D, E) up 
to 300 K and beyond (previously published data for the 
TEP of LMO [28] fall between those of specimens A2 
and C). This behavior and the smaller Nernst coefficient 
for specimen A2 are consistent with a departure from 
electron-hole symmetry in the more oxygenated crystals 
and incomplete cancellation of very large partial ther- 
mopowers. 

A two-component analysis of the transport places these 

qualitative observations on a quantitative footing. Mag- 

netoresistance and Hall data [18] for crystal C yield 
R h ^ _ R e ^ 3 4 x lcr s n m / T at 3 00 K and a h _ a e 

such that the measured Hall coefficient (Fig. 2) is an 
order of magnitude smaller than R h or |i? e |, consis- 
tent with the substantial cancellation implied by the 
thermoelectric coefficients. The effective carrier density 
ps 2/{R h \e\) ~ 3.7 x 10 26 nT 3 represents - 13% of that 
estimated from the chemistry/bonding [18], suggesting 
localization of a significant fraction of carriers at the FS 
in accord with optical studies [12]. A similar analysis 
of transport data for crystal A2 [18] yields R h ~ 7|i? e |, 
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a h ~ 0.4cr e , and a mobile carrier fraction ~ 2.5%. Us- 
ing these hole and electron transport coefficients and the 
measured v c b(T), the equation for v [25] can be inverted 
to compute S£(T) and S%(T) for the two crystals (inset, 
Fig. 3). In generating these curves it was assumed that 
cr' l /(T c = cr£/crb, but relaxing this assumption yields com- 
parable TEP magnitudes [18]. The colossal TEP values 
(2 — 6 mV/K at low T) are typical of semiconductors [29], 
but are unprecedented for metals. 

We hypothesize that the computed TEP is substan- 
tially larger than the intrinsic chain- axis value because 
the measured Ob in LMO is suppressed by the occurrence 
of incomplete electrical connectivity of the Mo-0 chains. 
This hypothesis is motivated by the Hall data implying 
a large fraction of localized carriers, and the observa- 
tion of a large constant term in «&(T) (12.9 W/mK) at 
T > 80 K [inset, Fig. 1(d)]. The latter is attributed to 
electronic heat conduction. Poor connectivity will sub- 
stantially suppress the measured Ob from that of contin- 
uous chains, whereas electronic heat conduction on chain 
fragments will continue to contribute given that heat is 
transferred readily to the lattice. Supporting this pic- 
ture, the constant term in Kb varies little among sev- 
eral other crystals with varying Gb [18]. This scenario is 
appealing because it can simultaneously account for the 
large fraction of total charge that is localized and the dif- 
ferent values for the anisotropy ratio ert>/er c reported in 
the literature [6, 16, 17, 30]. It also favors a Sommerfeld 
value for the Lorenz number (L = 2.45 x 10~ 8 Wtt/K 2 ): 
the resistivity for specimen C and its estimated mobile 
carrier fraction (/ = 0.13) from the Hall measurement 
yields, 4(300 K) = (L /f) x 0.55 x lO-^fhn)" 1 x 
300 K=10 W/mK. A bipolar thermal conductivity [31], 
expected in a metal with overlapping bands at Ep, should 
also contribute < 1 W/mK. Bipolar heat flow should 
dominate the thermal Hall conductivity (Righi-Leduc ef- 
fect) of a closely compensated system, offering an alter- 
native explanation for recent observations of an appar- 
ent Wiedemann-Franz law violation [30] in LMO based 
on such measurements. An intrinsic chain resistivity 
at 300 K, p c h = fpb — 70 /ifkm and Hall mobility 
fx c h ~ (pb/.f) — 2.3 x 10~ 2 T _1 are inferred. The com- 
puted partial thermopowers (Fig. 3), depending inversely 
on <7b [25], should be reduced accordingly by the factor 
/• 

On the other hand, the measured thermoelectric coef- 
ficients reflect appropriate bulk averages, weighing volt- 
age contributions from short and long chain fragments 
according to local temperature differences. Thus the 
intrinsic chain thermoelectric behavior must be closest 
to that of the most conducting specimen C, having the 
largest value of v and smallest of the TEP (the most 
closely compensated). The deviation from particle- hole 
symmetry and increased carrier localization character- 
izing more oxygenated crystals presumably reflect the 
combined effects of charge doping and disorder associ- 



ated with oxygen defects [15]. The electronic system of 
LMO is most one-dimensional at high T where the ther- 
mal energy exceeds the transverse hopping energy (tj_) 
that determines the warping of the FS sheets [32]. This 
is the regime where a conventional phonon-drag mech- 
anism is least likely to explain the large thermoelectric 
coefficients since anharmonic phonon-phonon scattering 
should limit momentum transfer to the electron system, 
rendering phonon drag negligible. However, it is possi- 
ble the low dimensionality of the lattice [13] suppresses 
phonon-phonon scattering along the chains [33], thereby 
enhancing phonon drag in this regime. It is difficult to as- 
sess the possible relevance of Luttinger physics. Theoret- 
ical treatments of the thermopower for single Luttinger 
chains with impurity scattering [34] indicate a linear tem- 
perature dependence, inconsistent with the present re- 
sults. More realistic models including electron-electron 
Umklapp scattering and coupled chains (essential for the 
Nernst effect) have not been treated to our knowledge. 
At T < Tm-, where higher-dimensional behavior is mani- 
fested in unconventional density-wave order [12-14] and 
the appearance of superconductivity, large phonon-drag 
effects arc plausible, particularly near T = 20 K where 
both Kb and v exhibit maxima (Fig. 1). There the phonon 
mean free path A is ~ 10 /im, as inferred from kinetic the- 
ory using the measured specific heat [12] and an acoustic 
phonon velocity v ~ 3 km/s. This implies a phonon- 
drag TEP [35], S g ~ vA/fi e - ph T w 750 /LtV/K using 
Me-p/i = 2T _1 (estimated as 5 times the average fi/f 
for specimen C), in good agreement with the computed 
partial TEP value (corrected by the factor /). 

In summary, Lio.9MogOi7 is found to have a bipolar 
Nernst coefficient that is among the largest known for 
metals and far exceeds the expected contribution from 
carrier diffusion over a broad temperature range. Though 
phonon drag appears capable of accounting for this dis- 
crepancy at the lowest T, low-dimensional physics of the 
electronic and/or lattice systems may be important to 
a complete understanding of its enhanced thermoelectric 
coefficients. 
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Supplementary Material 



J. L. Cohn, B. D. White, C. A. M. dos Santos, and J. J. Ncumeier 



Properties of the crystals studied 



TABLE I: Lattice constants, oxygen content as determined by energy-dispersive x-ray spectroscopy, and Hall coefficient (at 
250 K) for the six Lio.gMo60i7 crystals studied in this paper. 



Specimen a (A) b (A) c (A) Wt.% O R H (10~ 9 fi m/T) 

Al 12.751(2) 5.521(2) 9.475(5) 32.1(3) 

A2 a 12.752(2) 5.520(2) 9.482(4) 3.3 

B 12.751(2) 5.520(2) 9.486(5) 30.6(4) 

C 12.753(2) 5.520(2) 9.491(4) 2.1 

D 12.753(2) 5.519(2) 9.488(4) 

E 12.751(2) 5.520(2) 9.489(3) 29.7(3) 



"crystal A2 is the same as specimen Al, but measured 18 months later with fresh electrical contacts. The implication is that crystals 
tend to lose oxygen over time. 

Transport measurements 



For thermoelectric measurements specimens were suspended in vacuum with one end attached to a Cu cold sink 
(silver epoxy) and the other affixed with a small metal- him chip heater (varnish). Longitudinal and transverse volt- 
ages were measured with Au leads attached with silver epoxy. A 25 /im-diameter chromel- const ant an differential 
thermocouple monitored the temperature difference across the sample and a separate chromel-constantan thermocou- 
ple monitored the average sample temperature (relative to the calibrated thermometer mounted in the cold sink). 
A high-precision, differential method was employed at hxed magnetic fields to determine the Nernst coefficient: at 
each field the heater current was cycled through three different values and voltages measured after a steady state 
was achieved. This process was repeated for reversed magnetic field. The Nernst signal at each held was determined 
from the slope of resulting linear Vh~^T x curves, with Vr = \V X y{B) — V xy (—B)]/2 the held-odd transverse voltage 
(Fig. S2 a). The Nernst signal was linear in magnetic field at all temperatures for B < 3 T (Fig. S2 b). Nonlinearity 
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at higher fields was evident at temperatures T < 35 K. The Nernst coefficient was computed from the low-field slope 
in these plots. The field dependence of the thermocouple [1] for temperatures T > 1 K is < 5%. 

The Hall coefficient was measured during the same experimental runs. Linearity of the Hall voltage in applied 
current and magnetic field was observed for T > 50 K. Some nonlincarity in field was observed at lower temperatures. 
The Hall coefficients of both specimens A2 and C were re-measured in subsequent experiments under isothermal 
conditions (varnished flat to the heat sink); good agreement was found with both sets of measurements. 

For thermal-conductivity measurements, heat losses associated with radiation and conduction through the leads 
were assessed in separate experiments with specimens suspended from their leads in vacuum. The data presented in 
Fig. 1 d have been corrected accordingly. These corrections amounted to 30% near 300 K for k c and 5% for Kb- At 
T < 100 K the corresponding corrections were < 6% and < 1%, respectively. 



V H =[V xy (+B)-V(-B)]/2 



N=dV H /d(AT x ) 
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Fig. S 2: Variation of the Nernst voltage with longitudinal temperature difference at fixed fields (a) and field dependence of 
the Nernst signal at fixed temperatures (b). 



The heat current in highly anisotropic crystals always has a small component along the nominally transverse 
direction during Nernst measurements. This induces a field-odd Nernst voltage across the longitudinal voltage leads. 
Thus, for an experiment with heat flow along b, v\, c can be determined by extracting both the field-odd and field-even 
parts of V xx and AT y , respectively, the latter determined with an additional transverse thermocouple. This procedure 
yielded excellent agreement between values for v c b extracted from separate experiments on specimen C (discussed in 
Fig. 1) with heat currents J q \\c and Jq||6. A field-odd component of AT y produced by the Righi-Leduc effect was 
small, typically < 10 _2 AT K , consistent with the small Hall angle. 



Two-carrier analysis of magnetotransport 

The magnetoresistance (MR) and Hall resistivity of a conductor with two carrier types and field-independent 
densities and mobilities can be expressed as [2, 3], 



—— and PH (B) = — — 

1 + (PeffB) 2 1 + {HeffB) 2 



where, 



_{a h R h -a e R e ) 2 _a h a e {R h + R e ) 

MRo ° ~ a h o-e(R h + R e ) 2 ' Mi ~ o^T7 e 



{a h ) 2 R h + (o- e ) 2 R e R h R e 
Rho = , h , — j Rhoo — 



At T > 50 K the MR was quadratic and the Hall resistivity linear in field (Fig. S3), thus Rhoo is irrelevant. 
Simultaneously fitting the data at each temperature, with the additional condition that the measured zero-field 
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Fig. S 3: Magnetoresistance and Hall resistivity versus magnetic field for current applied along the conducting chains (b axis) 
and field along a. 
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Fig. S 4: Hall coefficients and conductivities for holes and electrons along the conducting chains (6 axis) from two-band model 
fitting. Error bars at T < 50 K reflect uncertainties associated with deviations of the MR from the two-band form (noted 
previously in Ref. 14). 

conductivity be equal to a h + <r e , highly constrains the four transport parameters. The reliability of parameter sets 
(determined using the Levenberg-Marquardt method) was established by varying the initial values of the parameters 
over a broad range. The results for crystals A2 and C are shown in Fig. S4. 

Fermi temperature and mobilities 

The Fermi temperature can be estimated from the energy difference between the Fermi level and bottoms of the 
two bands dispersing along the T-Y direction of the Brillouin zone, as determined by photoemission (0.25 and 0.49 
eV, Ref. 8) or band structure calculations (0.38 eV and 0.62 eV, Ref. 9), i.e. T F ~ (4 - 6) x 10 3 K. The dash-dot-dot 
curve labeled "283^/Tp" in Fig. 1 a uses Tp = 5 x 10 3 K and \i = /i' 1 + fi e computed as /J? = a l R l (i = h,e) from 
Fig. S4. 
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Fig. S5: Low-T v /T plotted versus [j./Tf for specimen C and a variety of other compounds (adapted from Ref. 5). 



Estimate of Fermi temperature from the thermopower 



Changes in the electronic structure at T < T m lead to a divergence of the TEP (Fig. 3) and thus preclude an estimate 
of Tp from the T — >• behavior of the diffusion thermopower. Instead we estimate Tp from the slopes of the linear-T 
thermopowers for specimens C, D, and E, under the assumption they reflect the diffusion contribution (we ignore the 
positive offsets in the TEP extrapolated to T — for these specimens that are reminiscent of the TEP in cuprates 
[4]). The average value, dS/dT = -0.045//V/K , combined with the Mott expression under the assumption of a linear 
energy dependence of the conductivity at the Fermi level yields, Tp = (7r 2 /3)(fcs/e)(— dS/dT) -1 ~ 6.3 x 10 3 K. This 
agrees satisfactorily with the Tp estimates above given that the linear slope of our data should be smaller than those 
of either electrons or holes due to partial cancellation, S = (a h S h + a e S e )/(a + a e ). 



Carrier density estimated from chemistry /bonding 

The two xy bands at Ep in the unit cell of Lio.9MoeOi7 (composed of two formula units) are half full (Ref. 9), i.e. 
1 electron per band or 2 mobile electrons per unit cell. The Li deficiency depletes charge (0.2 electrons per cell) from 
the chain Mo atoms, only half of which contribute to occupied orbitals at Ep, Thus the mobile charge per unit cell 
is estimated as 2 - 0.2/2 = 1.9 electrons. With the unit cell volume, V = 668 A 3 , this yields n = 2.84 x 10 27 m~ 3 . 
Oxygen deficiency should increase n from this value. 



Seebeck coefficients computed from anisotropic Nernst expression 

As noted in the text (Ref. 27), the Nernst coefficient for an anisotropic, bipolar conductor can be written as, 



(SI) 

where R h (R e ) and S h (S e ) are the hole (electron) Hall and partial Seebeck coefficients, respectively, and <r x = a^ + cr%, 
etc. for a y . The Seebeck coefficient is defined as, S x = (cr x S x + a x S x )/a x . For perfect compensation, <j x = a x , 
°y = <Ty, R h = —R e , the second term reduces to, (l/2)n x (S x — S x ). Note that a x = a x docs not necessarily imply 
dy = dy. Since we have not measured the magnetoresistance for electric current along the c axis, for the calculations 
of S£ and S% using eq. (SI) described in the text, we took <j^/<j c = l{< J bl (J b) with 7 as an adjustable parameter. In 
addition to the case 7 = 1 (shown in the inset of Fig. 3) we also show in Fig. S6 cases where holes are dominant for 
c-axis conduction (7 = 1.8, 3 for specimens A2, C) or electrons are dominant (7 = 0.1). Colossal Seebeck coefficients 
arc found in cither case. 
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Fig. S 6: Hole and electron partial thermopowers computed using equation (SI) and the two-band thermopower expression, 
with the measured Nernst coefficients and two-band fit parameters from Fig. S4 as input. Error bars (25%) are propagated 
from Fig. S4. 
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